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^ ! The acceleration of the expansion of the universe, and the need for Dark Energy, were in- 
^ ferred from the observations of Type la supernovae (SNe la)^"^. There is consensus that SNe 

I ; la are thermonuclear explosions that destroy carbon-oxygen white dwarf stars that accrete 

■ matter from a companion star^, although the nature of this companion remains uncertain. 
c« ■ SNe la are thought to be reliable distance indicators because they have a standard amount 
• • ! of fuel and a uniform trigger — they are predicted to explode when the mass of the white 

. ^ ; dwarf nears the Chandrasekhar mass" — 1.4 solar masses. Here we show that the high red- 
^ ; shift supernova SNLS-03D3bb has an exceptionally high luminosity and low kinetic energy 

■ that both imply a SMper-Chandrasekhar mass progenitor. Super- Chandrasekhar mass SNe 
la should preferentially occur in a young stellar population, so this may provide an explana- 
tion for the observed trend that overluminous SNe la only occur in young environments-'-. 
Since this supernova does not obey the relations that allow them to be calibrated as standard 
candles, and since no counterparts have been found at low redshift, future cosmology studies 
will have to consider contamination from such events. 

SNLS-03D3bb (SN 2003fg) was discovered on April 24, 2003 (UT) as part of the Supernova 
Legacy Survey (SNLS). Its redshift is 2; = 0.2440 ± 0.0003, determined from host galaxy [Oil], 
[OIII], Ha, and H/5 emission lines. A finding chart and observational details can be found in the 
supplementary information (SI). From the lightcurve (Figure [B we measure a peak magnitude in 
the rest frame V band, V = 20.50 ± 0.06 mag. This corresponds to an absolute magnitude of 
Mv = -19.94 ± 0.06 (Ho = 70 km s^^ Mpc-\ = 0.3, flat universe). SNLS-03D3bb faUs 



1 



completely outside of the My distribution of low-z SNe la-, and is 0.87 mag (a factor of 2.2) 
brighter than the median. Note that neither changes in the Hubble constant nor ^Im significantly 
affect this brightness difference. Asphericity may account for variations in SN la luminosity at the 
25% level, but not a factor of two^-*^. SNLS-03D3bb also does not follow the lightcurve width- 
luminosity relationship '° for SNe la that allows them to be calibrated as standard candles — it is 
too bright for its lightcurve width ("stretch", s=1.13) by 0.61 ± 0.14 mag (4.4o-). 

Type la supemovae are powered exclusively by the decay of ^^Ni and its decay product 
56Coii, requiring ~ 0.6 Mq of ^'^Ni to reproduce a normal SNe laimiii^. Since SNLS-03D3bb 
is 2.2 times overluminous, this implies that it has ~ 1.3 MQof ^^Ni. Such a large ^^Ni mass is not 
possible if the progenitor is limited to the Chandrasekhar mass (1.4 Mq). Even models that burn 
the entire 1.4 M© to nuclear statistical equilibrium via a pure detonation produce only 0.92 M© 
of ^^Ni, with the remainder comprising other iron-peak elements'^. Since at least 40% of the SN 
la must be elements other than ^^Ni to reproduce observed spectra^^^iii, this implies a WD mass 
of ~ 2.1 Mq. Some authors find rapid rotation may support such a massive white dwarf—. The 
merger of two massive white dwarfs could also produce a super-Chandrasekhar product-^^i^S. 

This simple estimation of the nickel mass is supported by a more detailed calculation us- 
ing the principle that the luminosity at maximum light is proportional to the instantaneous rate 
of radioactive decay^^^^. The impied Ni mass is^^: Mmi = , , where Lkoi is the bolo- 

ab{tii) 

metric luminosity at maximum light (the luminosity integrated from the ultraviolet to the in- 
frared), and a is the ratio of bolometric to radioactivity luminosities, near unity. S is the ra- 
dioactivity luminosity per solar mass of ^^Ni from its decay to ^^Co and subsequent decay to ^^Fe: 
S = 6.31 X lO'^^g-tfl/s.s _^ ;l_43 X iQ'is^-tn/iii g-i j^-i^ ^^Yere Ir is the time in days for 

the supernova to rise from explosion to maximum light. Using tji = s x 19.5 days^^, for SNLS- 
03D3bb, tji = 22 days (see SI for the effect of a shorter rise). We use a = 1.2 as a conservative 
value, although for high ^^Ni masses, a may be lower, since nickel above the photosphere will not 
contribute to the luminosity—. 

To convert our V magnitude into a bolometric equivalent, we use a synthetic spectrum cal- 
culated to match the observed UV-i-optical spectrum (Fig.|3l), but extended into the infrared (13% 
of the bolometric luminosity is from the IR extrapolation). The bolometric correction (mbc) is 
0.07 ± 0.03 mag, such that Mboi = My + = —19.87 ± 0.06 mag. Using these numbers, we 
calculate Mni = 1.29 ± 0.07 M© for SNLS-03D3bb, in agreement with the simple scaling argu- 
ment used earlier. The quoted error is from the statistical, k-correction, and bolometric correction 
errors added in quadrature. SNLS-03D3bb has a significantly larger bolometric luminosity and 
implied ^^Ni mass compared to low redshift SNe (Fig.EJ- 

SNLS-03D3bb also has an unusually low ejecta velocity, as shown in the Keck spetrum taken 
2 days after maximum light (Fig.|31l. With a Sill velocity of 8000 ±500 km s~^, it falls well outside 
the range of velocities seen for this feature at maximum light (Fig. HI). This is hard to understand 
in the Chandrasekhar mass model, which predicts higher velocities for more luminous SNe la, in 
contrast to the unusually low velocities in SNLS-03D3bb. 

The kinetic energy {E^) of a SN la arises from the difference between the nuclear energy 
(En) obtained from the synthesis of elements via fusion in the explosion and the binding energy 
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(Eh) of the white dwarf'^. Thus the kinetic energy velocity is: fke = y^M^^' where Mwd is 
the mass of the white dwarf. The binding energy of a 1.4 Mq C/0 WD is 0.5 x 10^^ ergs^"^. For a 
2M0 WD and a central density of 4 x 10^ g cm~'^, the binding energy is 1.3 x 10^^ ergs 

Since there are only three classes of elements in a SN la (iron-peak elements, intermediate 
mass elements (IME), and unbumed carbon and oxygen), a simple model can be developed for 
the nuclear energy generation. En. Burning a mixture of equal parts carbon and oxygen to the 
iron peak produces Epc = 1-55 x 10^^ erg Mq^^ while the synthesis of ^^Si produces 76% as 
much energy". Thus: En = -EpeAfwDl/Fe + 0.76/ime), where Mwd is in solar masses, and /pe 
and /iME are the fractional compositions of iron peak and intermediate mass elements. If /pc and 
/imp do not sum to one, the remainder is the fraction of unbumed carbon and oxygen (/c), which 
does not contribute to the nuclear energy. The ^^Ni makes up approximately 70% of iron-peak 
elements so we adopt Mni = 0.7Mwd/fc, where is the mass of ^^Ni. 

In the Chandrasekhar mass model, more luminous SNe, with more ^^Ni, have a higher t>ke 
(Fig.|4|). Increasing the fraction of unbumed carbon and oxygen, fc, can lower the kinetic energy, 
perhaps accounting for some of the dispersion in SN la velocities. However, this also lowers the 
available ^^Ni, so it cannot account for the low velocity seen in SNLS-03D3bb. 

The kinetic energy gives the velocity of the supernova averaged over the entire mass, ap- 
proximately equivalent to the velocity a few weeks after maximum light". The most appropriate 
observational signature of this velocity is unclear, since SN la line velocities change with time, 
and different ions can have different relative velocities. We find good agreement between the 
Sill velocity at 40 days after maximum light and the theoretical kinetic energy velocity, but we 
emphasize that this is in imperfect comparison. 

A super-Chandrasekhar mass reproduces the low velocities seen in SNLS-03D3bb (Fig.|4|). 
Since Chandrasekhar models with more Ni produce higher velocities, the low velocities of SNLS- 
03D3bb imply an increased progenitor binding energy and thus a larger total mass. As a caveat, 
we note that this simple calculation is only intended to illustrate general trends. Future theoretical 
studies will have to assess such complications as using different ions, different white dwarf density 
structures, and a wider range of binding energies. 

Super-Chandrasekhar mass SNe la should be more likely in a young stellar population, where 
the most massive stars exist-i^. The low mass, star forming host of SNLS-03D3bb is consistent 
with this scenario (see SI). Thus, the apparent existence of super-Chandrasekhar mass SNe la may 
explain why the most luminous SNe la only occur in young stellar environments^. The standard 
Chandrasekhar mass model offers no explanation for this behaviour, since the total amount of fuel 
and triggering mechanism should be independent of the mass of the progenitor stars. 

SNe such as SNLS-03D3bb will have to be screened out in cosmological studies. Since 
younger stellar environments produce more luminous SNe, as the mean stellar age decreases with 
redshift the mean properties of SNe la will change^. This can be calibrated if all SNe obey the 
same stretch-luminosity relationship, but SNLS-03D3bb does not. Its peculiarity was so obvious 
that it was excluded from the SNLS cosmological result-, but less extreme objects could lurk in SN 
samples. Future cosmology studies will have to carefully scmtinise SNe la from young populations 
to see if they obey the same lightcurve shape-luminosity relationship as other SNe la. 
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Figure 1: The lightcurve of SNLS-03D3bb. We fit k-corrected template lightcurves to tlie ob- 
served pliotometry of SNLS-03D3bb, tlien transform tlie peak magnitudes back to the Johnson- 
Cousins^ BV magnitudes in the Vega system. We find peak magnitudes of B = 20.35 mag and 
V = 20.50 ± 0.06 mag from a simultaneous fit to g' and r' data. The error (s.d.) consists of 0.04 
statistical error and 0.04 k-correction error. A lightcurve template was fit using the stretch method i 
(stretching the time axis of a template lightcurve by a stretch factor, s = 1.13). The epoch of max- 
imum light relative to the rest frame B band was determined from a simultaneous fit to all of the 
data. At maximum, we only use the V band value to compare to other SNe, since it is the best 
constrained. Data past -1-35 days was not used in the fit. The arrows are three sigma upper limits. 
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Figure 2: Bolometric luminosity vs. implied Ni mass for SNLS-03D3bb and low redshift SNe 
la^. The low redshift SNe la were fit using the same techniques as those used for SNLS-03D3bb: 
the bolometric luminosity was determined using the peak magnitude in the V band from a simul- 
taneous fit to B and V band data. For the low redshift SNe la, we integrated the s=l SN template 
to obtain a bolometric correction of 0.06 magnitudes and an uncertainty (s.d.) of 0.05 mag for 
the combined bolometric and k-correction error. The solid line represents a normal stretch=l SN 
la, with a rise time (U) of -19.5 days, while dotted lines show s = 0.9 (U = 17.6) and s = 1.1 
(tr = 21.5). Low luminosity SNe la have lower stretches, and thus shorter rise times, resulting in 
less ^^Ni for a given luminosity, while high stretch SNe la show opposite behaviour. The dashed 
line shows an upper limit for the expected ^^Ni mass in a Chandrasekhar mass SN la, obtained by 
burning the entire white dwarf to iron-peak elements in a detonation—. 
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Figure 3: Keck LRIS spectrum of SNLS-03D3bb at 2 days after maximum light compared to a 
spectrum of the normal Type la SN 1994D. Also plotted is a SYNOW fit to the data with dom- 
inant ions labeled. SYNOW is a parameterized resonance-scattering code, allowing the user to 
adjust optical depths, temperatures, and velocities to aid in the identification of supernova lines 
SYNOW parameters are listed in the supplementary information. SNLS-03D3bb shows the lines 
of IMEs typically seen in a SN la at maximum light — Sill, SII, and Call, but in SNLS-03D3bb 
the velocity of the lines is lower than usual. The line at 415 nm appears to be CII, but the other 
predicted carbon features cannot be clearly identified due to the lower signal-to-noise ratio of the 
spectrum in the red. No other identification could be found for the 415 nm feature. 
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Figure 4: Inferred Ni mass vs. Sill 615 nm velocity, a, The data at maximum light °. Ni masses 
are converted as described in the text using Mboi = Mb + 0.2. Red circles are from early type 
galaxies (E or SO), while blue circles are from late type galaxies (Sa-Irr). Circle size is proportional 
to fsi/^^si, where vsi is the rate of change of the velocity of the Sill feature. There is no vsi 
measurement for SNLS-03D3bb. b, Kinetic energy velocity of SNe la versus Ni mass for 1.4 
solar mass models with different fractions of unburned carbon (/c). This unburned fraction should 
not be much higher than ~ 20% because carbon is rarely seen in SN la spectra^". Overplotted 
symbols are vsi for low redshift SNe la^^ extrapolated to 40 days after maximum (correcting for 
stretch). For SNLS-03D3bb we use vsi from its closest neighbor. The error bar reflects the range 
if an average value of vsi is used. SNLS-03D3bb is not consistent with the 1 .4 Mq model, c, As 
above, but showing that A/wd ~ 2 Mq models can explain SNLS-03D3bb. Less extreme super- 
Chandrasekhar mass models are consistent with the low redshift data. The three low ^^Ni SNe are 
not necessarily super-Chandrasekhar SNe la — their large values of vsi make projections to 40 
days uncertain. o 



